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We present new lattice results of the K — ► -kIv semileptonic form factors obtained from simulations 
with two flavors of dynamical twisted-mass fermions, using pion masses as light as 260 MeV. Our 
main result is /+(0) = 0.9560(84), which, combined with the latest experimental data for decays, 
leads to \V U s\ = 0.2267(5)exp. (20) /+(0) . Using the PDG(2008) determinations of \V ud \ and \V ub \ our 
result implies for the unitarity relation \V ud \ 2 + \V U s\ 2 + \V u b\ 2 = 1.0004(15). For the 0(p 6 ) term of 
the chiral expansion of /+(0) we get A/ = / + (0) - 1 - f 2 = -0.0214(84). 

PACS numbers: 11.15.Ha,12.15.Hh,12.38.Gc 



The Cabibbo's angle, or equivalently the CKM matrix thus reduced to the problem of finding a prediction for 



element |V^ S | is one of the fundamental parameters 
of the Standard Model. The most precise determina- 
tion of \V US \ comes from K — > irtv (Kn) decay. The 
PDG(2008) quotes = 0.2255(19) % It is based 

on the new, very accurate experimental determination of 
the product |K s |/+(0) = 0.21668(45) @, H and on the 
old estimate of the vector form factor at zero-momentum 
transfer /+(0) = 0.961(8) given in Ref. Q. 

The determination of /+(0) using lattice QCD started 
only recently with the quenched calculation of Ref. 
where it was shown how /+(0) can be determined at 
the physical point with a ~ 1% accuracy. The findings 
of Ref. [|| triggered various unquenched calculations of 
/+(0), namely those of Refs. 0, S H with Nf = 2 and 
pion masses above ~ 500 MeV and the recent one of 
Ref. Q with Nf = 2+1 and pion masses starting from 
330 MeV. 

In this Letter we present a new lattice result for /+(0) 
obtained from simulations with two flavors of dynamical 
twisted-mass quarks, using pion masses from 260 MeV up 
to 580 MeV. Our determination of /+(0) includes the es- 
timates of all sources of systematic errors: discretization, 
finite size effects (FSE's), g 2 -dependence, chiral extrapo- 
lation and the effects of quenching the strange quark. 

The chiral extrapolation and the related uncertainty 
are investigated using both SU(3) and, for the first time, 
SU(2) Chiral Perturbation Theory (ChPT). Within the 
former one can perform a systematic expansion of /+ (0) 
of the type ,/ + (0) = 1 + / 2 + fa + where /„ = 
7r /(47r/ 7r ) Tl ] and the first term is equal to unity due 
to the current conservation in the SU(3)_ limit. Because 
of the Ademollo-Gatto (AG) theorem [lfj] , the first cor- 
rection /2 does not receive contributions from the local 
operators of the effective theory and can be computed un- 
ambiguously in terms of the kaon and pion masses {Mk 
and M„) and the pion decay constant It takes the 
value f% = —0.0226 at the physical point [4|. The task is 



A/ = / 4 + / 6 + ... = / + (0)-(l + / 2 ) 



(1) 



Recently SU(2) ChPT at the next-to-leading order 
(NLO) has been applied to study the quark-mass depen- 
dence of /+(0) [11]. In SU(2) ChPT the strange quark 
field does not satisfy chiral symmetry and the depen- 
dence on the strange quark mass, m s , is absorbed into 
the low-energy constants (LEC's) of the effective theory. 
The convergence of SU(2) ChPT is expected to be good 
when the u/d quark mass is significantly smaller than 
m s . In the case of /+(0) one gets the NLO result [llj 

q M 2 M 2 
/+(0) =F+- I (4^M^) + c+M 2 + 0{M$) (2) 

where F+ and c+ are LEC's functions of m s and c+ de- 
pends also on the renormalization scale [i in such a way 
that the whole NLO result $2$ is independent on /i. 

For the extrapolation of our lattice data to the physi- 
cal point we apply both SU(2) and SU(3) ChPT obtain- 
ing consistent results, which help constraining the uncer- 
tainty of the chiral extrapolation. 

We perform simulations with Nf = 2 flavors of dy- 
namical twisted-mass quarks [l2| generated with the 
tree-level Symanzik improved gauge action at a lat- 
tice spacing a = 0.0883(6) fm [H, Q (P = 3.9), 
for six values of the (bare) sea quark mass, namely 
am sea = 0.0030,0.0040,0.0064,0.0085,0.0100,0.0150 
(see Ref. [lj|). The valence light-quark mass is always 
kept equal to the sea quark mass (unitary pions) and the 
simulated pion masses goes from ~ 260 to ~ 575 MeV. 
For each pion mass we use three values of the (bare) 
strange quark mass, namely am s = 0.015,0.022,0.027, 
to allow for a smooth, local interpolation of our results 
to the physical strange quark mass (am p s hys ~ 0.021). 

At the two lowest pion masses the lattice volume is 
i 3 • T = 32 3 • 64 in lattice units, while at the higher ones 
it is 24 3 • 48 in order to guarantee that M-^L > 3.7. 
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We perform two additional simulations: the first one 
at M n ~ 300 MeV using the smaller volume and the 
second at M„ ~ 470 MeV using a finer lattice spacing 
a ~ 0.07 fm ((3 = 4.05) in order to check FSE's and 
discretization errors, respectively. 

The 2- and 3-point correlation functions relevant in 
this work are calculated using all-to-all quark propaga- 
tors evaluated with the "one-end-trick" stochastic proce- 
dure. All the necessary formulae can be easily inferred 
from Ref . [3] , where the degenerate case of the pion form 
factor is illustrated in details. At each value of the pion 
mass the statistical errors are evaluated with the jack- 
knife procedure, while a bootstrap sampling is applied in 
order to combine the jackknives for different pion masses. 

The matrix element of the weak vector current can 
be written as 

{n(p')\V,\K(p)) = P, f + (q 2 ) + q, f-(q 2 ) , (3) 

and the scalar 



where = p^ H 
form factor fo(q 2 



P't, and % = 
is defined as 



~Pu 



fo(q 2 ) = f+(q 2 ) + 



Ml 



f-(q 2 



(4) 



Following Ref. 
^2 



the scalar form factor at q 2 — q^ ax — 
(M K — M^) 2 can be calculated on the lattice with very 
high statistical precision using a suitable double ratio of 
3-point correlation functions. In the present simulations 
we get a precision better than ~ 0.2% (see Table Q]). 

At each pion and kaon masses we determine both the 
vector f+(q 2 ) and the scalar fo(q 2 ) form factors for sev- 
eral values of q 2 < q^ax m or der to interpolate at q 2 = 0. 
We take advantage of the twisted boundary conditions 
13j for details) to achieve values of q 2 quite 
= 0. The momentum dependencies of both 



(see Ref. 
close to q 



form factors are nicely fitted either by a pole behavior 

/+,o(g 2 )-/ + (0)/(i- s+ , o g 2 ) (5) 

or by a quadratic dependence on q 2 

f+Al 2 ) = /+(0) ■ (1 + S+,o q 2 + C+.o q 4 ) , (6) 

where the condition /o(0) = /+(0) is understood. The 
quality of the two fits is illustrated in Fig. [1] 

The values obtained for /+(0) depend on both the 
pion and kaon masses. The dependence on the latter 
is shown in Fig. [2] at ~ 435 MeV and it appears to 
be quite smooth, so that an interpolation at the phys- 
ical strange quark mass can be easily performed using 
quadratic splines. This is obtained by fixing the combi- 
nation (2M|- — M 2 ) at its physical value, which at each 
pion mass defines a reference kaon mass M^ 3 ■ 



2[M r K ef } 2 - M 2 = 2[M p K hys } 2 
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Fig. 1: Scalar fo(q 2 ) and vector f + (q 2 ) form factors obtained 
at Mir ~ 300 MeV and M K ~ 580 MeV versus q 2 in physical 
units. The solid and dashed lines are the results of the fits 
based on Eqs. and respectively. 
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Fig. 2: Results for /+(0) versus M 2 K at M w ~ 435 MeV. The 
square corresponds to the value of /+ (0) obtained by local in- 
terpolation via quadratic splines (dotted line) at the reference 
575 MeV from Eq. 0. 



kaon mass M 7 ^ 



Note that at the SU(3)-symmetric point Mr = M v 
the absolute normalization /+(0) = 1 is imposed auto- 
matically by the double ratio method of Ref. [B| . 

The results for /+(0), obtained using the pole dom- 
inance © or the quadratic fit ©, and interpolated at 
the reference kaon mass ([7]) , are given in Table [1] for each 
pion mass. It can be seen that the values of /+(0) corre- 
sponding to different ^-dependencies of the form factors 
differ by less than half of the statistical errors. In what 
follows we will show in the figures only the results ob- 
tained using the pole fit ((5|). 



(MeV) 


M r K et 
(MeV) 


/o(<?max) 


/+(0) 
(pole) 


/+(0) 
(quadratic) 


260 


520 


1.03097 (224) 


0.97519 (499) 


0.97374 (505) 


300 


530 


1.01923 (121) 


0.98052 (440) 


0.97950 (390) 


375 


555 


1.00961 (123) 


0.98916 (264) 


0.98813 (248) 


435 


575 


1.00416 (43) 


0.99343 (130) 


0.99273 (131) 


470 


590 


1.00272 (34) 


0.99421 (79) 


0.99413 (85) 


575 


635 


1.00016 ( 6) 


0.99823 (15) 


0.99827 (19) 
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Table 1: Results for fo(qmax) and f+(0), obtained with the 
pole 0) or quadratic ^Bj) fits, interpolated at the reference kaon 
mass £?P f° r each simulated pion mass. 



with Ml h v s = 135.0 MeV and M 



phys 
K 



494.4 MeV. 



The SU(3) chiral analysis of /+(0) starts by considering 
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the NLO term using the exact expression evalu- 
ated for our partially quenched (PQ) setup in Ref. 

2M 2 K + Ml ZM 2 K Ml\og(M 2 jM 2 K ) 



f PQ 

J2 



32tt 2 j 'I 6Air 2 fZ(M 2 K -M%) 
M\{±Ml - M2)log(2 - Ml/Ml) 



(8) 



6i**ft(M% -Ml) 

and by constructing the quantity A/ from Eq. (fT]). We 
then carry out the extrapolation to the physical point 
using a simple phenomenological ansatz in terms of Ml : 

Af = A + AtMl + A 2 M* + A 3 Ml log(M 2 ) , (9) 

where Ao^i^s are fitting parameters. 

The results obtained for /+(0) using two fits for A/, 
one with A3 = and the other with A2 = 0, are shown 
in Fig. [21a). It can be seen that: i) the (absolute) size 
of A/, whose chiral expansion starts from the NNLO 
term /<4, is even larger than the one of the leading NLO 
term / 2 at all pion masses, and ii) the impact of the 
logarithmic term at NNLO is quite small. 
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Fig. 3: Results for /+(0) versus Ml at Mk = M^ analyzed 
using SU(3) (a) and SU(2) (b) ChPT. In (a) the SU(3) LO 
+ NLO term, l + f^®, is shown by the dashed line. The solid 
and dot-dashed lines are the results of the fit for Af with 
A3 = and A2 = 0, respectively. In (b) the dashed line is 
the SU(2) LO + NLO fit ^ applied to our data with < 
0.4 GeV, while the solid line corresponds to the result of fitting 
all lattice points adding to Eq. a NNLO term proportional 
to Ml. The vertical line corresponds to M% hys = 135.0 MeV. 



A relevant check on our fits © is that they turn out to 
be consistent with zero (within the statistical errors) at 
the point M^ — M r ^ , as required by the AG theorem. 

At the physical point we get 



/+(°)&) = 0.9599(61)(32) 



where the first error is statistical and the second one 
is systematic coming from the uncertainties of the mass 
extrapolation and the g 2 -dependence of the form factors. 

We now discuss the analysis based on SU(2) ChPT. 
First we note that Eq. © holds for full QCD as well 
as for the PQ theory with Nf — 2. In the latter case it 
can be verified by expanding f£® [see Eq. ([8])] in powers 
oiMl/M^. Thus we consider a SU(2) fit of the form © 
treating F + and c + as fitting parameters, and we apply it 
to our data with M T < 0.4 GeV. Alternatively we add to 
Eq. @ a NNLO correction proportional to Ml and apply 
the new fit to all lattice points. The results are shown in 
Fig. CJb). It can be seen that the impact of the SU(2) 
NNLO correction is quite small up to M w w 0.5 GeV at 
variance with the corresponding SU(3) result shown in 
Fig. a). This finding signals a better convergence of 
SU(2) ChPT with respect to SU(3) for /+(0). 

At the physical point we get 



/+(°)ls(7(2) 



0.9563(53)(13) 



(11) 



(10) 



The application of SU(2) and SU(3) ChPT yields re- 
sults for /+(0), Eqs. (fTT7|) and (fTl"]) . which are consistent 
within the uncertainties. By averaging the two results 
and adding the systematic errors in quadrature we get 

/f 3 (0) - 0.9581 ± 0.0057 stat . ± 0.0035 syst . . (12) 

We now present our estimates of the remaining sources 
of systematic effects. 

Finite Size. We have performed a simulation at M n = 
300 MeV using the volume 24 3 • 48, which corresponds 
to M W L ~ 3.2. We get /+(0) = 0.98633(362) using the 
pole-dominance fit ([5]) and /+(0) = 0.98597(337) using 
the quadratic fit ([6]). We combine these values with the 
results shown in the second row of Table [TJ corresponding 
to the volume 32 3 • 64 with M^L ~ 4.2. Assuming a 
volume dependence of the form A + Be~ Ml < L / L 3 / 2 we 
obtain a residual FSE equal to 0.0018, which we add (in 
quadrature) to the systematic error of Eq. ([12"]) . 

Discretization. We have performed a simulation at 
Afjr ~ 470 MeV using a finer lattice spacing (a ~ 
0.07 fm). We observe a systematic increase of the 
scalar form factor fo(q 2 ) at all values of q 2 and for all 
kaon masses. In particular we get /+(0) = 0.99555(80) 
and 0.99518(95) using the pole-dominance ([5]) and the 
quadratic © fits, respectively. We combine these val- 
ues with the results shown in the fifth row of Table [1] 
Assuming a linear fit in a 2 (which is consistent with the 
automatic 0(a) improvement at maximal twist [13])) we 
find a discretization effect equal to 0.0037, which we add 
both to the central value and (in quadrature) to the sys- 
tematic error of Eq. (fl2|) . Clearly a more detailed study 
of the scaling property of /+(0) would be beneficial in 
order to estimate better and to reduce further the dis- 
cretization error. 



4 



Quenching of the strange quark. The effect of our PQ 
setup can be estimated within SU(3) ChPT, because, 
thanks to the AG theorem, the effect of quenching the 
strange quark is exactly known at NLO: at the physical 
point h - f[ Q = -0.0058 (~ 26% of f 2 ). This correc- 



tion is added to the central value of Eq. (|12[) . As for the 
0(p 6 ) term Af, we have found evidence that the chiral 
logs, which are the most sensitive to quenching effects, 
are small compared to the contribution of the local terms 
(see Fig.[3](a)). We estimate that the quenching effect on 
Af is at most 50% of the same effect on f 2 . Thus we 
add (in quadrature) the value 0.0028 (~ 13% of Af) to 
the systematic error of Eq. (TT2")) . Note that this value is 
of the same size of the difference between our estimate of 
Af at Nt = 2 and the quenched one of Ref . [5J . 
Our final result is 



/+(0) = 0.9560 ±0.0057 stat .± 0.0062 syst . 
= 0.9560 ±0.0084 , 



(13) 



which corresponds to Af = —0.0214(84). Our determi- 
nation agrees very well with the Leutwyler-Roos result 
[3| and with previous lattice calculations at Nf = 
N f = 2 0,0,11 and N f = 2 + 1 @|. 

Using the latest experimental determination of the 
product |K s |/+(0) = 0.21668(45) HQ we get from p]) 



|K, | = 0.2267 ±0.0005 exp . ± 0.0020 /+(0) 



(14) 



Combining this value with \V ud \ = 0.97418(27) and 
1Kb | = 0.00393(36) from PDG2008 [2] the CKM uni- 
tarity relation becomes 



IKdl 2 



IKJ 2 



IK,! 2 



1.0004 ±0.0015. 



(15) 



In conclusion we present our results for the slopes of 
the scalar (so) and vector (s + ) form factors. Their light- 
quark mass dependence is illustrated in Fig. [4] and it 
appears to be quite mild. We have tried simple fitting 
functions of the form 



djM* log(M*) 



(16) 



where a,j, bj, 



and dj are fitting parameters and j — 
+, 0. The results of two fits, one with d +i o = and the 
other with c +i n = 0, are shown in Fig. |4] 

In terms of the dimensionless quantities A + .n = 
M% s +i o the extrapolation to the physical point and the 
evaluation of the systematic uncertainties yield 



An = 



(12.8±2.2 stat .±4.5 syst .)-10- 3 , 
A+ = (23.7 ± 2.3 stat . ± 2.1 syst .) • 10~ 3 , (17) 

where the large systematic error on An is dominated by 
discretization effects. Our results for both An and A+ 
agree very well with the latest experimental averages 
X e xp - = (13.4 ± 1.2) • 10- 3 and A+ p ' = (24.9 ± 1.1) • 10~ 3 , 
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Fig. 4: Results for the slopes so (dots) and s+ (squares) versus 



M* at M K = M 



ref 



The solid (dashed) line corresponds to 



the fit S16V with d+,o = (c+,o — 0). 



obtained in Ref. using data from KLOE, KTeV, IS- 
TRA+ and NA48 experiments. 

We thank all the ETMC members for fruitful discus- 
sions and the apeNEXT computer centres in Rome and 
Zeuthen for their invaluable technical help. One of us 
(F.M.) also acknowledges the Consolider-Ingenio 2010 
Program CPAN (CSD2007-00042). 



[9. 
[10 

[11 

[12. 

[is; 

[14 

[is; 
[is; 

[17 



N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531; M. 
Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 
(1973) 652. 

C. Amsler et al. [Particle Data Group], Phys. Lett. B 667 

(2008) 1. 

M. Antonelli et al. [FlaviaNet Working Group on Kaon 

Decays] , larXiv:0801.T8T7l [hep-ph] . 

H. Leutwyler and M. Roos, Z. Phys. C 25 (1984) 91. 

D. Becirevic et al, Nucl. Phys. B 705 (2005) 339. 

N. Tsutsui et al. [JLQCD collaboration], PoS LAT2005 
(2006) 357. 

C. Dawson, T. Izubuchi, T. Kaneko, S. Sasaki and 
A. Soni, Phys. Rev. D 74 (2006) 114502. 

D. Brommel et al. [The QCDSF collaboration], PoS 
LAT2007 (2007) 364. 

P. A. Boyle et al, Phys. Rev. Lett. 100 (2008) 141601. 
M. Ademollo and R. Gatto, Phys. Rev. Lett. 13 (1964) 
264. 

J. M. Flynn and C. T. Sachrajda [RBC/ UKQCD collab- 
oration], Nucl. Phys. B 812 (2009) 64. 
R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz [Alpha 
collaboration], JHEP 0108 (2001) 058. 
R. Frezzotti, V. Lubicz and S. Simula, Phys. Rev. D 79 

(2009) 074506; 

P. Dimopoulos et al, [ETM collaboration], 
larXiv:0810.2873l [hep-lat]. 

Ph. Boucaud et al. [ETM collaboration], Phys. Lett. B 
650 (2007) 304; Comput. Phys. Commun. 179 (2008) 
695. 

D. Becirevic, G. Martinelli and G. Villadoro, Phys. Lett. 
B 633 (2006) 84. 

R. Frezzotti and G. C. Rossi, JHEP 0408 (2004) 007. 



